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1 Introduction 



The fractal structures of random geometrical objects have been under intensive 
investigation for a number of years, both in connection with quantum gravity pQ 
and in the study of disordered materials [21131 E| This work is to a large extent aimed 
at understanding the geometric characteristics of generic objects in the ensembles 
under study and how these characteristics are reflected in physical phenomena. 

An important notion in the study of fractal geometries is the concept of dimen- 
sion. Definitions of dimension which agree on smooth manifolds do in general not 
do so in the random or fractal case. One important concept is that of spectral di- 
mension which is defined to be d s provided the heat kernel at coinciding points, 
averaged over the random geometries and viewed as a function of time t, decreases 
as t~ ds / 2 as t — > oo. Equivalently, the spectral dimension is a measure of how likely 
a random walker is to be at the starting point after time t. This notion of dimension 
is in general different from that of Hausdorff dimension djj which is defined in terms 
of the growth of the expectation value of the volume of a geodesic ball of radius r 
as r — > oo: 



We will study in detail many examples of this for the case of random combs in this 
paper. The discrepancy between these dimensions is also well demonstrated, at least 
numerically, in quantum gravity. 

Early work on the spectral dimension in quantum gravity was done by numerical 
simulation [Sj. This lead to the investigation of random walks on random trees 
and the spectral dimension of random trees was calculated analytically in the 
extension to non-generic trees was given in In [H] a scaling relation was derived 
which relates the spectral dimension to the extrinsic Hausdorff dimension. Related 
work on the spectral dimension of trees can be found in EI]- in the condensed 
matter community the spectral dimension has been investigated for a variety of 
systems, see, e.g., j2]. While very few exact results have been obtained, random 
combs in particular have been studied numerically as well as by mean field theory 
methods. Mean field theory simplifies the problem since it allows one to model the 
walk on a comb by a random walk on the spine of a comb with a waiting time 
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distribution which is taken to be the same for all vertices on the spine. The waiting 
times arise from the excursions that the walk makes into the teeth of the comb. The 
spectral dimension of random combs with teeth whose lengths obey a power law 
distribution has been studied in [H], see also ^2]- If the exponent of the power law 
distribution is a then the spectral dimension was found to be given by 



for a < 2 but 1 otherwise. In this paper we prove this result which shows that mean 
field theory is exact in this case. Mean field theory was shown to be exact in a 
special case in We will also show that the spectral dimension of random combs 
whose teeth may be infinitely long is always 3/2. 

We develop technical tools to prove the above mentioned results and also apply 
those techniques to random trees and to some examples of non-random combs. The 
main new idea is a splitting of random walks as well as random combs into subsets 
that either yield exponentially suppressed or uniformly controllable contributions to 
the quantities under consideration (typically d s ). The tools are recursion relations 
for generating functions, simple monotonicity results and convexity arguments. We 
believe that the methods can be applied to study random walks on more complicated 
random graphs. 

In the next section we define the random comb ensembles we wish to study and 
the most important generating functions and critical exponents. We establish simple 
monotonicity results and use them to obtain some elementary bounds. In Section 
3 we study random combs which have an infinitely long tooth at each site on the 
spine with a nonzero probability. In this case the spectral dimension is always 3/2 
which is the same as the spectral dimension for a comb with all teeth infinitely long. 
In Section 4 we calculate the spectral dimension of combs with random but finitely 
long teeth and show that the spectral dimension is determined by the tail of the 
length distribution. In Section 5 we apply our methods to prove upper and lower 
bounds on the spectral dimension of random trees. In Section 6 we calculate the 
spectral dimension of fixed combs whose toothlength increases along the spine and 
also combs with infinite teeth whose separation increases along the spine. Section 
7 contains results about transport along the backbone of the combs and the full 
heat kernel on random combs. In the final section we discuss the relevance of our 
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methods and results for random geometry in general and compare our work with 
relevant results in the mathematics literature. Various technical calculations are 
relegated to two appendices. 



2 Preliminaries 

Let iVoo denote the nonnegative integers regarded as a graph so that n has the 
neighbours nil except for which only has 1 as a neighbour. Let N# be the integers 
0, 1, . . . , t regarded as a graph so that each integer n G Nf has two neighbours nil 
except for and £ which only have one neighbour, 1 and £ — 1, respectively. A comb 
C is an infinite rooted tree-graph with a special subgraph S called the spine which 
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Figure 1: A comb. 



is isomorphic to with the root at 0. At each vertex of S, except the root 0, 
there may be attached one of the graphs N% or N^. We adopt the convention that 
these linear graphs which are glued to the spine are attached at their endpoint 0. 
The linear graphs attached to the spine are called the teeth of the comb, see Fig. 
1. We will find it convenient to say that a vertex on the spine with no tooth has a 
tooth of length 0. We will denote by T n the tooth attached to the vertex n on S, 
and by Ck the comb obtained by removing the links (0, 1), . . . , (k — 1, k), the teeth 
Ti , . . . , T k and relabelling the remaining vertices on the spine in the obvious way. 
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2.1 Random walks on combs 

We consider simple random walks on the combs. We assume that the walker starts at 
the root unless we specify otherwise. At each time step the walker steps with equal 
probabilities to one of the neighbouring vertices. This means that the walker has 1, 
2 or at most 3 choices of vertices to step to at any given time and the corresponding 
probabilities are 1, 1/2 and 1/3. We are interested in the asymptotic properties of 
the walk after many time steps. We regard time as integer valued. 

Given a comb C, let pc{t) be the probability that the walker is at the root at 
time t. Let Pc(t) be the probability that the walker is at the root for the first time 
after t — at time t with the convention Pc(0) = 0. Clearly both pc{t) and p]j{t) 
vanish unless t is an even number. Given a random walk u which comes back to the 
root at time t it is clear that this may be the first return, the second one, etc. We 
can therefore decompose pc{t) into a sum over walks that have had a fixed number 
of intermediate visits to the root before ending there at time t, i.e., 

oo n 

Pc(t)=5 t , + J2 E UpcM- ( 3 ) 

n=l t 1 +t 2 +...+t n =t j=l 

We define the generating functions for return to the root and first return to the root 

by 

oo 



Qc(z) = >Jpc(t) (4) 



t=0 



and 



It follows then from (|2J) that 



P c (z) = J2^Vc(t). (5) 



The function Pc{ z ) is clearly analytic in the unit disc and satisfies |P<7(z)| < 1 for 
\z\ < 1 . Note that functions analogous to Pc and Qc are defined for the simple 
random walk on any rooted graph V with root of order one. We will denote these 
by Pr and Qr- If T consists of a single vertex, i.e., V = N Q , we adopt the convention 
that Pv{z) = 1. We will see in a moment that -Pc(l) = 1 f° r ah combs, i.e., the 
simple random walk is recurrent on combs. 
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Consider a fixed comb C. Any walk that contributes to Pq can be decomposed 
into a first step from to 1, then an arbitrarily large number of round trips into the 
tooth Ti intermingled with round trips into the comb C\ and then a final step from 
1 to 0. Each time the walk is located at 1 the probability of stepping into T\ or C\ 
is 1/3 and likewise the probability of the final step to is 1/3. It follows that 

p °<*> = (?) 

In particular, if C has no tooth at 1, we have 

Pc{z) 



z 2 



In fact the recurrence relations (J7J) and (jHJ) are valid for generalized combs where 
the teeth are allowed to be arbitrary rooted graphs (with a root of order 1), not 
necessarily the linear ones that we study here. 

Consider now the toothless comb N^. The generating function for first return 
to the root, denoted P^, satisfies 

It is convenient to introduce the variable x related to z by 

l-x = z 2 . (10) 

The generating functions are even functions of z so they can be regarded as functions 
of x; we will denote them by the same symbol which should not cause confusion and 
assume from now on that < x < 1. From Q we see that 



Pooix) = 1~VX. (11) 

For a finite tooth Ng we denote the generating function for first return of random 
walks to by Pg. An elementary calculation using the recurrence relation (jHj) and 
P\{x) — 1 — x yields 

see Appendix 1. We observe that Pe{x) is a decreasing function of i for a fixed x 
and PAx) — > Poo(x) as £ — > oo. 
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The comb which has an infinite tooth at each vertex on the spine will be called 
the full comb and denoted *. By (J7J) the function satisfies 

P * (X) = 3-pX)-P.(x)' (13) 
Using (|TT|) and the fact that P*(x) < 1 we find that 



P # (x) = 1 - + ^+1^. (14) 

We define the critical exponent a for a comb C by 

1 - P c (x) ~ x Q , as x -> 0, (15) 

where f(x) ~ ^(x) means that for any e > there are positive constants c% and C2 
such that 

C\X e fix) < g(x) < c 2 x~ £ f(x) (16) 

for < x < 1. We see from (jllj) and (fT4*j) that the half line and the full comb have 
a = 1/2 and a = 1/4, respectively. It is easy to relate a to the spectral dimension 
d s . If 

Pc (t)~r d °/ 2 (17) 

as t — ► oo then 

Qofc) ~ a;" 1+ds/2 (18) 

as x — »• and 

4 = 2- 2a, (19) 

so the half line and the full comb have spectral dimensions 1 and 3/2, respectively. 
The value d s = 3/2 for the full comb was first obtained in jHj. In fact the spectral 
dimension (if it exists) of any comb lies in the closed interval [1,3/2]. This is a 
consequence of 

P*(x) < P c (x) < P^x) (20) 

valid for any comb C . The inequalities fl2()D follow from the Monotonicity Lemma 
below. Furthermore, the lower bound in (|2()jl and (|14j) imply that random walks on 
combs are recurrent as claimed above. 

We note from (JZj) that, for fixed x, Pc(x) is a monotonic increasing function of 
-PtiOe) and Pc x {x). By applying (7) in turn to Pc\ (x), Pc 2 (x), . . .Pc k _ L (x) we find 
by induction the following result: 



Lemma A The function Pc(x) is a monotonia increasing function of Pt 1 (x), . . . 
P Tk (x),P Ck (x) for any k>l. 

Monotonicity Lemma Pc(x) is a decreasing function of the length, i^, of the tooth 
Tj, for any k > 1 . 

Proof By Lemma A Pc{x) is a monotonic increasing function of Px k {x) which is a 
decreasing function of according to (JT2J) . 

Rearrangement Lemma Let C be the comb obtained from C by swapping the 
teeth T n and T n+ \. Then Pc{x) > Pc(x) if and only if PT n {x) > Pr n+1 (x). 

Proof By Lemma A, and noting that T 1; . . .T n _i are the same for C and C, it 
suffices to prove that -Pc n _i > Pc' _ x if an d only if Pt„ > Pr n+1 - By (J2J) the former 
holds if and only if Px n + Pc n > Pr n+1 + Pc< ■ It is therefore enough to compute 

A = P Tn + Pc n - Pr n+1 ~ Pc n 

= (Pr n ~ Pr n+1 ) ( 1 - p p p p r ) (21) 

where we have used (JJJ) and the fact that C n+ i = C' n+1 . We see that A > if and 
only if Pr n > Pt„ +1 which completes the proof. 

2.2 The two-point function 

Let C be a comb and let Pq(P, n) denote the probability that a random walk that 
starts at the root at time is at the vertex n on the spine at time t and has not 
visited the root in the time interval from to t. We will refer to the generating 
function for these probabilities as the two-point function and denote it by Gc(x; n). 
Note that G c (x; 0) = P c (x) and 

oo 

G c (x;n) = ^(l-a;f^; n ). (22) 
t=\ 

The two-point function can also be expressed as the sum over all random walks from 
to n which avoid and end at n: 

M-i 

G c (x;n)= ]T JJ (a{uj t )- l yfT^;) (23) 

u:0-*n, ujt^O if t^O t=0 

where \u\ is the number of steps in the walk u and o~{uj t ) is the order of the vertex 
uj t where u is located at time t. 
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The representation (|2Hj) of the two-point function is quite useful and allows us 
to relate it to the first return generating functions as we now show. Let C be a 
comb and let Ck be defined as before. If we consider a random walk wonC which 
contributes to the two-point function Gc{x] n) we can decompose it into a sequence 
of n random walks uj 1 , . . . , uj n where uj 1 is a walk from to 1 which is identical to uj 
until uj leaves the vertex 1 for the last time before going to n, uj 2 is a walk from 1 to 
2 which is identical to uj after it left 1 for the last time until it leaves 2 for the last 
time, etc. The last walk uj n is the part of uj after it left n — 1 for the last time and 
until it ends at n. If for each k = 1, 2, . . . , n — 1 we add a last step to uj k back to 
the vertex k — 1 and call the resulting walk uj k we see that uj k is a walk from k — 1 
to k — 1 which contributes to Pc k _ 1 and any walk contributing to -Pc fc _i can arise in 
this way. It follows that for n > 

n-1 

G c (x; n) = a(n)(l - x)^ 2 J] P 0k (x), (24) 

fc=0 

where a(n) is the degree of the vertex n. 

We see from from and (J2U that the two-point function for the half line, 

G OQ (x;n), is given by 

nil 

X x 



G.(x ! »)=2^ rj -^J (25) 

for n > 0. We define the mass, m(x), of the two-point function Gc by its rate of 
exponential decay, i.e., 

logG c (a:;n) , v 

mix) = — hm . (26) 

For an arbitrary comb there is no reason to expect the limit (|26J) to exist but the 
mass associated with the two-point function for the half line is clearly 

7^ (a?) - -log 1 + V Z . (27) 



We can similarly use (|T4j) to compute the two-point function, G*(x; n) and the mass, 
m*(x), for the full comb. It furthermore follows from (|2(Jj) and (J24|) that 

^'"^ G„(x;n) < G c (x;n) < (^-) G^fon) (28) 



for any comb C. If the mass m(x) exists we define its critical exponent v by 

m{x) ~ (29) 
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as x — > 0. It is easy to see that v = 1/2 for and v = 1/4 for to*. From (|28j) 
we conclude that the critical exponent of the mass for any comb lies in the interval 

U' 2 J" 

The above considerations show that the exponents a and v coincide for the half 
line and the full comb. Indeed, we will prove in Section 6 that the scaling relation 

a = v (30) 

holds quite generally. 

2.3 Random combs 

Let C denote the collection of all combs. Let denote the nonnegative integers. 
If we are given a probability measure \i on U {oo} we can define a probability 
measure ir on C by letting the length of the teeth be identically and independently 
distributed by \x. This means that the measure of the set of combs Q with teeth at 
rii, n 2 , • • • , n k having lengths £ h £ 2 , . . . , 4 is 

k 

7r(n)=IjM^)- (31) 

j=l 

We will refer to the set C equipped with the probability measure 7r as a random 
comb. Measurable subsets A of C are called events and ir(A) is the probability of 
the event A. We define the first return generating functions for random walks on 
random combs as 

P(x) = (P c (x)) (32) 
where (•) denotes expectation with respect to the measure ir, i.e., 

(F{C)) = J F{C) d7r (33) 

for any 7r-integrable function F defined on C. Similarly, 

Q(x) = (Qc{x)}. (34) 

If Q(x) ~ x~ 1+ds l 2 as x — > we say that the spectral dimension of the random comb 
is d s . Similarly we define the exponent a for the random comb by 1 — P(x) ~ x a . 
We will see for the examples of random combs studied in this paper that the relation 
CH| holds. 
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The two-point function of the random comb is defined as 

G(x; n) = (Gc(x; n)). (35) 



We show below that the mass exists for any random comb. It follows from (|28j) that 
this mass rh(x) satisfies the inequalities 

moo(x) < fh{x) < m*(x) (36) 

and the critical exponent v of m lies in the interval [7, |]- 

2.4 The mass for random combs 

In this subsection we introduce some auxiliary generating functions and prove the 
existence of the mass for random combs. We assume that we are given a random 
comb where the lengths of the teeth are identically and independently distributed. 
For a fixed comb C define a modified two-point function Gq(x; n) by restricting the 
sum in f!23[) to walks that stop the first time they hit the vertex n. Then we have 
the factorization 

G c (x; n) = Gq(x; n)Q c (x; n), (37) 

where Qc{x;n) is the sum over all walks which begin and end at n and avoid the 
root 0. Equation (|37j) can be obtained by considering any walk contributing to the 
two-point function Gc{x;n) and cutting it at n the first time it hits n. The first 
part then contributes to Gq(x; n) while the second part contributes to Qc{x\ n)- Let 
T n be the tooth of C at n (which may be empty). Let Pq \x; n) be the generating 
function for first return of random walks that begin at n, have a first step to n — 1, 
and avoid the the root. Similarly let Pq(x\ n) be the generating function for first 
return of random walks that begin at n and have a first step to n + 1. Then we have 

Q c {x; n) = n . (38) 

3-Pt\x;n)-Pt\x;n)-P Tn (x) 

Using Pq\x] n) < 1 — y/x, Pq \x; n) < 1 and Pr n { x ) < 1 we obtain 

Qc{x;n) < 3x^5 (39) 

and hence, 

G° c (x; n) < G c (x; n) < 3x"3G^(x; n). (40) 
11 



Let us define G° c (x\ n, n') as the sum, analogous to ([23)1 . over all walks from n to n' 
which avoid both n and n' at all intermediate times. Then Gq(x; 0, n) = Gq(x; n). 
Consider now a walk u contributing to the two-point function Gq(x; ri\ + 712). Cut 
u the first time it hits n\. Cut it again the last time it leaves n\. Then we obtain 3 
walks, the first of which contributes to G° c (x\ rii), the second one starts and ends at 
7ii, avoiding both and rii + n 2 and the last one contributes to Gq(x; n\, rii + n 2 ). 
We therefore obtain a factorization 

G° c (x; ni + n 2 ) = G° c (x; ni)R c {x\ m)G° c {x\ ni, n x + n 2 ), (41) 

where Rc(x;n) < Qc{x;n). Hence, by (j39|) . 

G° c (x; ni + n 2 ) < 3x~^G° c (x; n x )G° c {x] n 1 ,n 1 + n 2 ). (42) 

Since the teeth are independently distributed we see that the functions G° c (x;ni) 
and G^(x;ni,ni + n 2 ) are also independently distributed. We denote the averaged 
modified two-point functions by G°(x;n,n'). Then 

G°{x;n 1 ,n 1 +n 2 ) =G°{x;n 2 ), (43) 

so 

G°(x; m + n 2 ) < 3x-^G°(x; n l )G°(x] n 2 ), (44) 

and the function log(3x~2 G°(x;n)) is subadditive in n. It follows by standard 
arguments that 

_ ^ \ogG°(x;n) = _ \ogG°(x;n) ^ 

n^oo n n n 

In view of (141])) we conclude that the mass associated with the averaged two-point 
function G(x] n) exists and is given by (JJ5J). 

2.5 Some bounds 

In this subsection we establish bounds which, for the purpose of calculating a and 
d s , allow us to ignore walks that wander too far along the spine. We denote by E(x) 
a generic nonnegative function of x > with the property that there are positive 
constants c%, c 2 and e such that 

E{x) < Cl e- C2X ~\ (46) 
12 



In the following we will let c, c', ci, C2 etc. denote positive constants whose value may 
change from line to line. 

Let C be a comb and define P^\x) as the contribution to Pc(x) coming from 
walks whose maximal distance from the root along the spine is n. Then 



oo 

,(n) 



P c (x) = J2 P c( x ) ( 47 ) 



n=l 



and we claim that 



if 



This follows from 



P^{x)=E{x) (48) 

n=N(x) 

G c (x;N(x)) = E{x). (49) 



Pc\x) = cr{N{x))- l GUx;N{x)fQ c {x ] N{x)) 

n=N(x) 

= a(N(x))- 1 G c (x;N(x)) 2 Q c (x;N(x))- 1 , (50) 

cf. (J3*7f) . (J3U|) and pPjl. We conclude, in particular, that if N(x) > x~^~ £ for some 
e > 0, then (HHJ holds for any C. 

Finally, if we have a comb C such that (J49j) holds for C and that C" is another 
comb which is identical to C up to the vertex N(x) on the spine, then, by ()37|) and 

(HOD, 

oo 

]T pW(i) < (G° c ,(x;N(x))) 2 Q c ,(x;N(x)) 

n=N(x) 

< cx^(G° c (x;N(x))) 2 

= E(x). (51) 
The estimates (jlH| and (jKTjl will be used repeatedly in this paper. 

3 Combs with infinite teeth at random location 

In this section we consider a random comb for which there is probability p G (0, 1) 
that there is an infinite tooth at a vertex on the spine and probability q = 1 — p that 
there is no tooth, i.e., /i(0) = p and n(oo) = 1 — p in the notation of subsection 2.3. 
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We will show that in this case the spectral dimension is 3/2, i.e., the same as for the 
full comb. It follows immediately that any random comb with a nonzero probability 
of an infinite tooth at any given vertex has spectral dimension 3/2. 

The strategy of the proof is to use (|4*H|) which shows that for a given value of x 
it suffices to consider walks that do not move beyond a location N(x) on the spine. 
Then we use the Rearrangement Lemma to dilute the teeth on the interval from 
to N(x) = [x~^~ £ ] so that they are regularly spaced and we can obtain an upper 
bound on Q which turns out to be of the same form as the trivial lower bound on 
Q coming from comparison with the full comb. 

Let Lq denote the distance from the root to the first (non-trivial) tooth and let 
Li, i > 1, denote the distance from the ith tooth to the (i + l)st tooth. Since the 
Lj's are independently distributed random variables we see that 

ndLi < L : i = 0, 1, . . . , k - 1}) = (1 - q L ) k . (52) 

If r is a real number we denote its integer part by [r]. Fix e G (0, 1/8), choose k = 
[x~"2~ £ ] and L = [x~ e ]. Let A e be the event that > L for some i G {0,1,..., k}. 
Then, by ir(A £ ) = E{x). 

Now consider a comb C ^ A £ . The spacings between the first k teeth of C are 
all smaller than or equal to L. By removing all teeth in C except the first k ones 
and shifting these suitably away from the root we obtain a comb C whose teeth 
have constant spacing Li = L, i = 0, . . . , k — 1. Hence, by the Monotonicity and 
Rearrangement Lemmas we have 

Pc(x) < Pc>{x) = P { c){x) + P®(x), (53) 

where Pq)(x) is the contribution to Pc'(x) coming from paths which do not pass 
through the point [a;~2 _£ ] on the spine and Pg) (x) is the remainder. By (}4*%|) we 
have Pq, (x) = E{x) uniformly for C G C \ A e . Moreover, we have 

pg,\x) = PS(x)< P* L (x), (54) 

where *l is the comb with infinite teeth of spacing L. We conclude from ()53)1 and 
dSH) that 

P c (x)<P* L (x) + E(x). (55) 
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Using the result 

P* L ( X ) < (56) 
derived in Appendix 2, we obtain 

Pc{x) < 1 - cx3+i (57) 

and it follows that 

Qc(x) < cx~^~^ (58) 

uniformly for C ^ A £ . Hence, 

Q(x) = [ Q c (x)dir(C)+ [ Q c (x)dn(C) 

J As JC\Ae 

< x~^n(A £ ) + cx~^~^n(C\A £ ) 

— 1 — £ 

< cx 4 2. (59) 

It follows that a < \ + | and d s > | — £ for any e > 0. In view of the lower bound 
(|2~U|) we obtain 

1 3 

a = 4' ds= 2' ^ 

4 Combs with finite random teeth 

In this section we will calculate the spectral dimension of random combs with finite 
but arbitrarily long teeth. An upper bound on P(x) will be obtained by mimicking 
the argument for the upper bound obtained in the previous section using the fact 
that if the teeth are sufficiently long they can be replaced by infinitely long teeth 
up to discrepancies of size E(x). The lower bound will be obtained by a convexity 
argument. 

Let \i be a probability measure on the non-negative integers and set /jl(£) = fie. 
For simplicity of presentation we assume to begin with that we have a power law 
distribution fi£ = c a £~ a , where a > 1 since the \ig sum to 1. Choose e > and define 

oo 

for x > 0. We shall refer to p(x) as the probability that a tooth is long. Clearly 
p(x) — > — > 0. More precisely, 

p{x) ~x 9{£ \ (62) 
15 



where 

g(e) = ^ + e(a-l). (63) 

Consider now the random comb defined by ji and let M = [x~ 9 ^~ e ] and K = 
[x~2~ £ ]. Denote by Lq the distance from the root to the first long tooth and by Lj, 
i > 1, the distance from the zth long tooth to the i + 1st long tooth. Let B £ be the 
event that at least one of the distances L , . . . , L^-i is greater than M. Since the 
probability that Lj > M is given by 

(l-p(x)) M <e- Mp{x \ (64) 

we have 

tt(B £ ) < Ke~ Mp{x) = E(x). (65) 

Consider now a comb C B e . By deleting all teeth from C except the first K 
long teeth and shifting these suitably away from the root we obtain a comb C whose 
teeth have a constant spacing Lj = M, i = 0, . . . , K — 1. By the Monotonicity and 
Rearrangement Lemmas we have 

Pc{x)<Pc>(x). (66) 

Replacing the teeth in C by infinitely long ones, only changes Pc'{x) by E(x). This 
follows from the fact that walks that go more than a distance [x~5~ e ] into the teeth 
only contribute E(x) to Pc{x) uniformly in C . By the same argument as the one 
leading to (jB3j) we obtain 

P c (x)<P* M +E(x), (67) 

uniformly for C ^ B £ . 

Let us first consider the case a < 2. Then we can choose e so small that g(e) +e < 
1/2. It follows that M ~ [x~ p \ with < (3 < 1/2 so by (12121) in Appendix 2, 

P c (i)<1-ctHW^) j ( 68 ) 

for C B e . We conclude from the above and (f63J) that 

< 1 - cx^ i90+£) (69) 
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for e and x sufficiently small, where we have defined go = g(0). Hence, 

«< i + y = i - (70) 

Similarly we conclude that 

Q(x) < cx~%. (71) 

If a > 2 the upper bound on a obtained above is replaced by the trivial upper bound 
a < | coming from the comparison with the comb with no teeth and similarly for 
the exponent of Q. 

We now turn to the proof of the lower bound. We first note that 

Q(x) > - \ (72) 
1 — Pyx) 

by Jensen's inequality. It will suffice to find a suitable lower bound on P(x). Noting 
that the lengths of the teeth are independently and identically distributed we take 
the expectation of (JJJ) to get 

Hx) = t»( *-rM-Poto )- (73) 

Applying Jensen's inequality again we obtain 

p (x) > o p\\ X p, v (74) 
3 — Pt(x) — P{x) 



where 



Pt(x) = J2^ P ?( x )- (75) 



1=0 



By rearranging the inequality (|74j) and using the fact that P(x) < 1 it follows that 



P(x)>l-yJl + x-Fr(x). (76) 
Using the expression (fT2"|) for Pe(x) and the definition (|27j) we can write 

_ l _ p—m{x)t 

1 - P t {x) = ^ 1 + e _ m(x)e (77) 
where we have put m{x) = 2m OQ (x) for convenience. It follows that 

oo 

1 - Prix) < X^Ml - e- em{x) ). (78) 

1=0 
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Consider first the case a > 2. Then, from (J7%|) . we obtain the bound 

oo 

1 - P T (x) < ^fx m(x) t-H < ex. (79) 

e=o 

Hence, P(x) > 1 — Cy/x and we conclude that a = | and d s = 1. 
More generally define 

oo 

I y = ^2n e P (80) 
e=o 

for 7 > 0. Let 

7o = sup{7 > : J 7 < oo}. (81) 

For the power law distribution \ig = c a £~ a we see that 70 = a — 1. For a < 2 it 
follows by a similar argument as before, using 

1 - e~ z < z p , z>0, < p < 1, (82) 

that for £ > 0, 

1 - PtO) < cv^m(x) 70 ~ e . (83) 

We conclude from (f7U|) that 

P(x) > 1 - ca;^^- 5 ) (84) 

so a > a/4. Combining this with ()70|) shows that a = a/4 and therefore the spectral 
dimension is given by 

ds = ^ (85) 

for 1 < a < 2. 

Let us now consider the general case when is an arbitrary probability distribu- 
tion and assume first that 70 defined as above is finite. For simplicity, let us further 
assume 5 that there exists a non-increasing function g(e) such that p(x) ~ x 9 ^ £ \ 
Then 

limg(e) = g = ^ . (86) 

eJ-0 Z 



5 This assumption is not true for arbitrary probability distributions and in that case the subse- 
quent arguments will have to be modified. We will leave this technical point aside. 
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In order to prove this let 7 > and consider the integral 



J — p(x)x~?' y ~ 1 dx 
Jo 
„1 00 

= / V HiX-^<- 1 

J0 !_ 



"2" 

00 r \ 



< 



/ , x-^dx. (87) 



Doing the x integral above we see that 



2 

J < -c + - /Je(£ + 1)1+^ (88) 

which is finite if 7 < (1 + 2e)7o. Similarly we can show that 

2 00 

J>-c (89) 

and the right hand side in the above inequality diverges if 7 > (1 + 2e)7 . We 
conclude that 

^) = ^7o(l + 2e) (90) 

so 9o — |to- This, together with the previous arguments, proves that d s — 1 if 
1 < 7o < 00 an d c? s = (3 — 7o)/2 if < 70 < 1. It is easy to check that d s = 1 for 
probability distributions which have 70 = 00. 

5 Random trees 

In this section we show how the results and methods of the previous sections can be 
used to bound the spectral dimension of infinite planar random trees. We begin by 
recalling some results about such trees from [THj . We consider rooted trees where 
the root has order 1. There is a measure 7r on these trees obtained as a limit of 
the uniform measures on ensembles of trees with a finite number of vertices. With 
respect to this measure there is with probability 1 a unique infinite simple path in a 
random tree r which can be viewed as a "spine" with finite trees attached to the 
spine. For simplicity of presentation, let us consider trees whose vertices have order 
1, 2 or 3. Then there is at most a single rooted tree Tj with root of order 1 attached 
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to each vertex j '• ^ on the spine of the infinite random tree. It is shown in ^Bj 
that the trees Tj are identically and independently distributed with the probability 
distribution 

p({rj}) = Z-^-N (91) 

where \tj\ denotes the number of links in Tj and Z is a normalization factor (which 
in this particular case happens to be equal to 1). 

The spectral dimension for the infinite random tree <i* ree is now defined in the 
same way as the spectral dimension for random combs by considering the return 
probability to the root for a simple random walk on the random tree. We will show 
that 

\ < dT < \. (92) 

The lower bound is obtained by showing that the generating function for first return 
of random walks to the root on the random tree, Pt Tee {x), is bounded from above by 
the corresponding generating function for random combs whose tooth length have 
a power law distribution with exponent 3/2. The upper bound is obtained by a 
convexity argument. The result (}9~2*|) is in agreement with d^ ee = 4/3 found in [S| 
by different methods. 

In order to prove ()92j) we first note that 

3 3 

Ci77T2 < p({ r ; | r | — 77,}) < c 2 n~2 ; (93) 

see, e.g., [Tj. By Lemma A and we see that it is sufficient for the lower bound 
in ([92)1 to show that among all the finite rooted trees r with a given number I of 
links, it is Ni which has the largest first return generating function. In order to 
prove this consider a tree r with I links and let P T (x) be its first return generating 
function. Let v be a vertex in r at a maximal distance from the root. Let uj be the 
unique simple path in r from the root to v. Then v has order one and r can be 
viewed as a comb with a finite spine uj and finite trees (possibly empty) attached 
to each vertex of uj. Let Vj, j = 0,1, ... ,n, be the vertices of uj, v n = v. Let tj 
be the tree attached to Vj. Let tk be the first non-empty tree we encounter as we 
move from v along uj towards the root. Let r' be the tree obtained by swapping tk 
and tk+i (which is empty by hypothesis). By the argument of the Rearrangement 
Lemma we deduce that P T (x) < P T i(x). Taking now a vertex v' in r' at a maximal 
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distance from the root and repeating the above argument we construct a sequence 
of trees and for some finite value i = i$ we obtain r^°) = N#. 

In order to prove the upper bound in ()92j) we begin by noting that by Jensen's 
inequality 

^tree(^) > 1 - , r (94) 

where 

P t {x) = ±523-Mp r {x) (95) 

T 

and the sum in (|95j) runs over all finite rooted trees with vertices of order 1, 2 or 3. 
By the same argument as the one leading to f)7fij) we find that 



Ptvee(x)>l- \Jl+X-P t (x). (96) 

We will now show that 

P t (x)>l-cy/x (97) 

which implies the desired result. Given a tree r it can always be decomposed into 
the root link and two subtrees ta and tb as indicated in Fig- El Note that ta and/or 




Figure 2: A decomposition of a tree into two subtrees. 
tb may be empty. This allows us to write 

/jx |ta|+|tb|+i l _ x 
Z 2 ^ I 3 J 3-P TA (x)-P rs (x)- (9? 



Using Jensen's inequality again and Z = 1 we obtain 

P t (a:) > -(1 -x) f 1 + 7tt^ + ^rr-r ) ■ (") 

3 V ' V 2-P t (x) 3-2P t (x)/ 
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Rearranging ((QSj) leads to 

(6P t (x) - ll)(P t (x) - l) 2 + x(2Pl(x) - WP t (x) + 11) > (100) 
which implies (19 7j) since Pt(x) — ► 1 as x — ► 0. 

6 Some non-random combs 

In this section we use the techniques we have developed so far to calculate the 
spectral dimension of some non-random combs. We first remark that it follows easily 
from the previous results that any periodic comb has either spectral dimension 3/2 
(if it has some infinite teeth) or spectral dimension 1 if all teeth have finite length. 

We will discuss two different types of combs: (i) combs with infinite teeth where 
the distance between tooth n and tooth n + 1 is an increasing function of n and 
(ii) combs with finite teeth such that the length of tooth n is an increasing function 
of n. For simplicity of presentation we will choose the increasing functions to be 
powers but our methods apply to more general cases. In both cases we find that the 
spectral dimension varies continuously with the power. 

6.1 Combs with increasing tooth spacings. 

Let a > and define C to be the comb all of whose teeth are infinite such that the 
distance from the fcth tooth to the k + 1st tooth is 

L k = [(k + 1)% (101) 

where [r] denotes the integer part of r as before. The distance from the root to the 
first tooth is Lo = 1. We will show that this comb has spectral dimension 

d s = (102) 
2 + a 

We first prove an upper bound on Pc(x). This will be done by an inductive 
argument bounding the two point function at a suitable distance along the spine: 

G c {x, [x-^]) = E(x) (103) 

for arbtrary e > 0, where 

a + 1 (104) 



2(a + 2) 
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For this purpose we need to introduce slightly more general combs than C . For any 
/ > 1 let C* be the comb with all teeth infinite such that the distance from the fcth 
tooth to the k + 1st tooth is 

L k =[(k + fn (105) 

Note that Lq = [f a ] in this case and C 1 = C. The inductive hypothesis is as follows. 
There exists r/ G (fj, |] such that for any r\ > r] and any constant c > there is a 
fixed E'-function, as defined in ()46|) . and the inequality 



G cf (x;[x-*>\) <E{x) (106) 

holds for all / in the range 1 < / < cx 1 + a . By (J49)) the hypothesis is true for 
T] = ~ since ()106j) holds in this case for any value of / > 1. We shall prove that the 
statement then holds also with r] replaced by 4>(r)o), where 

4>{rj) = -+ , ar] . . (107) 
yyu 4 2(o+l) K J 

The strategy of the following argument is to use the induction hypothesis ()106|) to 

obtain an upper bound on P C f(x) which in turn will be used to prove an improved 

upper bound on G c s using the representation (J2H). 

Consider one of the combs C* . The distance from the root to the nth tooth is 

given by 



k=0 

and fulfills 

_J_(( n + / _ i)-+i _ (/ _ i)«+i) < D n < 1 ((„ + - /•+!). (109) 

a + 1 a + 1 

For a fixed 77 > 770 choose n as a function of x > 0, such that -D n+ i > [x~ v ] > D n . It 

follows that 

L" 1 ~ x^+i (110) 

since / a+1 < c a+1 x _J? by assumption. The contribution to P C f(x) coming from 
walks that go beyond D n+i is bounded by E(x) as a consequence of (|106j) . By the 
Rearrangement Lemma and (|5ip it follows that 

P C f(x) < P* ln (x)+E(x) 

< l-cx v ', (111) 
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where rj = 4>(rj) and we have used (jllOJl and (j212j) to obtain the second inequality. 
In order to complete the inductive step we now need to show that ()106j) holds with 
7] replaced by rj . So let be given with 1 < / < cx~ v ' and let D = [x~ v ']. By 
and the Rearrangement Lemma we have 

D-l 



G cf (x,D) = tr(D)(l - x)-* D ft P c 

k=0 

< (l-x)^ D (P cf/ (x)) D , (112) 



where the comb is defined such that the distance from its root to the first tooth 
is larger than or equal to the largest tooth separation up a distance D along the 
spine in G* . For this purpose it suffices to take 



f = c 'x-^ (113) 

for a suitable constant c'. The comb is in the class covered by the induction 
hypothesis with rj replaced by r] — e for a suitable e since rf < 77. Hence, 

P cf >(x)<l- Cl x^- £ \ (114) 

and we conclude from (jll2|) that 

G cf (x,[x-*]) = E(x) (115) 

for 1 < / < cx _r? '. We have thus proven that if the induction hypothesis holds for 
a particular rj it also holds for 771 = 4>(r)o). Defining r\ r inductively by 

Vr+l = 4>{rir) (116) 

we see that {q r } is a decreasing sequence in the interval [ 2 (^. 2 ) ; 5] ari d the induction 
hypothesis holds for all 77 which satisfy 

n > fj = hm Vr = (117) 



r— »oo 



Pc(x) < 1 - cx^+^+i) + E(x) (118) 



We conclude that 
for any rj > fj and therefore 

JfcCaO < 1 - ciW +£ (119) 
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for any e > 0. This proves that 



< a + 1 ( } 

~ 2(a + 2) K J 



In order to prove the corresponding lower bound on Pc{x) let 77 G (77, i±^). Let 
_D n ~ [x~' n ] as before and let C be the comb C with all teeth beyond the rath one 
removed. Then from (j5Tj) and f)106j) it follows that 

P c ,(:r) = P c (x) + £(x). (121) 

Define C(n) to be the comb with infinite teeth at the vertices 1, 2, ...,n but no 
other teeth. Then, by the Rearrangement Lemma and (|121j) . 



o — • — • — • — • — • — • — • — • — • • • 

L 

Figure 3: The comb C(L) with L = 5. 



Pc(x)>P c(n) (x) + E(x). (122) 

The asymptotic behaviour of Pc(n)(%) as x — > is computed in Appendix 1 and we 
find 

Pc{x) > l-cxi~^T (123) 
which implies the desired converse to the inequality (|119|h 

Pc(x) > 1 - cx^~ e (124) 

for any e > 0. We conclude that 

a = " +1 (125) 

2(o + 2) v ; 

and the spectral dimension ()102|) follows. 
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6.2 Combs with increasing tooth length 

In this subsection we consider the comb C for which the length of is given by 
4 — [k a ], a > 0. We will prove that the spectral dimension of this comb is | if a > 2 
but 

2(1 + a) 

d s = v ; 126 
2 + a v ; 

for < a < 2. In order to prove this we first consider the case < a < 2 and define 
the comb C by 

4 = 0, k < k 

4 = [K], k>k , (127) 

where ko will be chosen to depend on x below. The comb C can be obtained from 
C by shortening teeth so by the Monotonicity Lemma 

Pc{x)<Pc>{x). (128) 

We then set ko = [x~P], < (3 < |, and choose j3 to optimize the bound ()128|1 . The 
asymptotic behaviour of the function Pc (x) is computed in Appendix 1 where this 
function is denoted Pk/{x) with k = fco and i = [k^]. 
Taking first the case /3a < \ we find 

P c ,[x) = 1 - ex 5 + 0(x s+£ ), (129) 

where 5 = max{ 1 ~ 2 ^ a , 0} and e > 0. The value of 5 is minimized by choosing 
P = Popt = 2+^ an d we conclude that 

« < 77^' ( 13 °) 

2 + a 

provided that /3 op ta < |; that is for a < 2. 

To obtain a lower bound on Pc{x) we first note that by monotonicity Pc k (x) < 
Pc{x), as Cf. can be obtained from C by lengthening teeth, and therefore 

n-l 

G c (x;n) = (1 - x)- n ' 2 a{n) J[ P Ck {x) 

k=0 

< 3(1 -x)~ n l 2 (P c (x)) n 

< c\ exp ( — C2 ni^ J . (131) 
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Combining this fact with (|HT|) and the Monotonicity Lemma shows that 



Pc{x) > P c »(x) +E(x), (132) 

where C" is the comb with teeth of constant length = [kf] and k\ = [x _5 +^~ £ ]. 
This comb has exactly the structure of the comb # considered in Appendix 1. We 
find, using p05jl . 

Pc(x) > 1 -cx^~ £ ' (133) 
for any e' > 0. It follows that a > combining this with (j!30p gives the results 

a = d s = + a , < a < 2. (134) 

2 + a 2 + a 

We now turn to the case a > 2. We use the argument leading to ()129j) with 
(3 — \. The teeth of C are now so long that they are effectively infinite and (| 193)1 . 
(gUDD and (EHTt yield 

^c(a:) < ^C'(x) = 1 -cx3 + 0(x3 +£ ) (135) 

where e > 0. From ()20|) we know that Pc^) — P*{x) and it follows immediately 
that 

a = - A , 4 = | «>2. (136) 

7 Anomalous diffusion 

In this section we explore the connection between the full heat kernel on random 
combs and the functions we have focused on so far, namely the two-point function 
and the first return generating function. The main result is that anomalous diffusion 
along the spine is described by the decay of the two point function and the critical 
exponents a and v coincide. We will focus on the random comb with random 
toothlength. The comb with random spacing between infinitely long teeth can be 
treated by similar arguments. 

7.1 The exponents a and v are equal 



Our starting point is the representation (j24j) of the two-point function on a comb 
C . We will bound v from below with a using a convexity argument. The opposite 
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inequality follows from pointwise estimates as obtained in Sections 3 and 4 in the 
calculation of the spectral dimension. 

Using (J2J) and (|24jl . and remarking that 

1 



> e~ y (137) 



l + y 

for any y > 0, we find 

n) > a(n)(l - x) n / 2 exp ( - £(2 - P Tfc+1 (a;) - Pc fc+ ») ) • (138) 

\ k=0 J 

Averaging over the comb ensemble and applying Jensen's inequality we obtain 

G(x; n) > (1 - x )n/2 e -n{2-P T (x)-P( X )) ( 13 g) 

where Pt{x) is the average of the first return generating functions on the individual 
teeth. Clearly Pt{x) > P(x) so 

G(X; U) > (1 - x )«/2 e -2n(l-P(x))_ ( 14Q ) 

We conclude immediately that v > a. 

In order to prove the converse inequality it is sufficient to show that 

G{x; [x- a - £ ']) = E(x) (141) 

for arbitrarily small e' > 0. Indeed, it follows from the definition of the mass and 
(JUJ) that 

G(x; n) > G\x- n) > e - m{x)n (142) 

for all n. Hence f!141|) implies m(x)x~ a ~ 6 ' — > oo as x — > 0, which shows that 
< a + e'. 

To establish (|141|) . we split the average over C into a contribution from B £ and 
a contribution from C\B e as we did in Section 4 for P(x). By ()28)1 the former 
contribution is bounded from above by 

^G 00 (x;n)7r(B £ ) = E(x). (143) 

For C £ C\B £ we recall from Section 4 that the first return generating function can 
be estimated by 

Pc(x) <P* M (x) + E(x), (144) 
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where M = [x~ 9 ^~ £ ], and E(x) is independent of C. We claim that the bound (jl44j) 
also holds for C^, k = 1, . . . , uniformly in C \ B £ for e' sufficiently small. 

Recalling that combs C G C\B £ are charcterized by the requirement Lq, . . . , Lk-i < 
M, where K = [x~z~ e ], and that 4 < a < |, we choose e' such that a + e' < | + e, 
and hence x~ a ~ s ' / K — > as a; — > 0. In particular, for each under consideration, 
we have that L 0} . . . , -^[i^j < M, for x sufficiently small, and inspection of the proof 
of (j!44|) shows that the inequality holds for the CVs as claimed (with a modified 
^-function) . We can thus use (jl44j) together with the product representation (}2~4"j) 
and obtain for the second contribution to G(x; [x~ a ~ £ ']) the bound 

(1 - x) n/2 (P* M (x) + E(x)) [x ~ a ~ £ ' ] = E(x) . (145) 

We have thus proved (jl41}) and hence also v < a. 

7.2 The heat kernel 

Consider a comb C and define Kc{t; n, k) as the probability that a random walker 
who leaves the root at time is located at the vertex k in the nth tooth of C at 
time t. We will denote this vertex as (n, k). If the nth tooth of C has length smaller 
than k we define this probability to be 0. We will refer to the function Ka(t] n, k) 
as the heat kernel on C since it satisfies the heat (or diffusion) equation on C: 

K{t + 1; n, k) = (T ( n '' K)- x K(t\ n, k'), (146) 

(n'.fc') 

where the sum in f!146|) runs over the nearest neighbours of the vertex (n, k) in C. 
Next define the function 

oo 

K c (t;n) = Y,K c (t;n,k) (147) 

fc=0 

which is the probability that a walker has travelled a distance n along the spine at 
time t. We are interested in the asymptotic behaviour of Kc{t; n) for large t and n. 
In order to analyse this function we define the corresponding Green function by the 
Laplace transformation 

oo 

H c (x; n) = ^(1 - xf 2 K c (t; n). (148) 
t=o 
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Decomposing the walks contributing to the heat kernel we can express Hc{x;n) in 
terms of previously defined generating functions as 

H c (x-n)= Gc{ ^ ] D e (x) (149) 
1 - P c {x) 

where i is the length of the nth tooth of C and 

1 1 

D e {x) = l + ~^2G Nl {x;k). (150) 

fc=i 

Using ({23| and (fT2"j) we can write Dg(x) in terms of Pg(x) as 

/ s 2 1 , , \ (\ - Pfi(x)\ 

D e (x) = - + -(l + VT^) [ — ^J. 



(151) 



Hence, 



We begin by establishing a lower bound on He- Let C(oon) denote the comb C 
with the nth tooth replaced by an infinite tooth. By the Monotonicity Lemma and 
it follows that 

3 

Pc(oon)(x) < P c (x), G C (oon)(x;n) < -G c (x;n). (152) 

(I _ x)~ n / 2 DAx) n_1 
H c (x;n) > i- 1 JJ Pc (a; ). (153) 

We note that Dp(x) is the only quantity on the right hand side of the inequality 
()15Hj) which depends on the length of the nth tooth of C. Now consider one of 
the combs Ck{oon — k) and swap teeth n — k — 1 times so that the infinite tooth 
becomes the first tooth while the tooth Tj of Ck{oon — k) becomes tooth number 
j + 1 for j = 1,2, . . . ,n — k — 1. Denote the resulting comb by C' k (oon — k). By the 
Rearrangement Lemma P Ck (oon-k)(x) > P C ' k (oon~k)(x). By (0), 

1 — X 

Pc' k (oon-k)(x) = ~ ~ n . (154) 

3-Poo(x) -P^+iW 

where C = C'(oon). We now average over C. Using (J137)) . combined with ()154j) . 
and Jensen's inequality we obtain, 

(H c (x; n)> > G-^<^ e -"^W-PW), (155) 

We now take the ensemble average of D^. Let us first assume that the teeth are 
finitely long with probability 1 and let 70 be given by (JSTj). Then, by (|S3*jl . 

1 - (P/(x)> < cx5( 1+ i»- £ ). (156) 
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The converse inequality, with e replaced by — e, follows from (jH§j) and (jSljl so 



1 - (P t (x)) ~ x 2c \ (157) 

We conclude from (|151|) that 

(D t (x)) > Cl + c 2 x 2a -\ (158) 

By (HH2D, (Pc(oon)(ac)) < On the other hand, (TH)3|) with k = and Jensen's 

inequality imply that 

S 3-P.W-P(x) 

and we conclude that 

l-P(ar)~l-<P c(oon) (x)). (160) 

Combining ffTKHl) . ifTSEjl and fWi|) yields 

(if c (x; n)) > ci /- 1 e- C3M ", (161) 

In order to establish the corresponding upper bound we use the inequalities ()144j) 
and ()145|) which imply that 

(I — x)~ n / 2 Di>(x) 

H c (x; n) < { } 1 ) ((P* M (x) + E 1 (x)) n + E 2 (x)(P 00 (x)) n ) , (162) 

-L — Pm{X) 

where the ^-functions £i and i?2 are uniform in C. It is clear from ()151j) that the 
converse of the inequality ()158|) holds (with different constants Ci and c 2 ). Hence, 

(ff c (z; n)) ~ f-V-^ 1 " + E(a;)e- n ^ (163) 

for small x. If infinite teeth appear with nonzero probability in the comb ensemble, 
then a = | and we can ignore the finite teeth since they do not affect the critical 
behaviour. In this case, ()163|) holds with a = \. 

The asymptotic relation (|163|) enables us to compute the mean extent down the 
spine of walks at large time. We have, for k > 0, 

"£(l-xf 2 ((n% :H=t ) c = 5>*<l?c(s;n)> 

t n 

+ E(x), x^O, (164) 
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and so, by standard Tauberian theorems (see, e.g., [201121]), 

( {n k )^\ =t )c ~ c 2 t ka , t -> oo. (165) 

On the other hand, ()163|) does of course not allow us to compute (K c (t; n)). How- 
ever, if we make the ansatz, 

(K c (t;n))^c^e- C '^, (166) 
we find that this is consistent with ()163|) if and only if 

1 — a 1 — a 2(1 — a) 2(1 — a) 

This is in agreement with the exact results for the half line [T2] and the full comb 
[IE]. We believe that ()166j) captures the essential behaviour of the averaged heat 
kernel for 

n 

— > 1. 168 

Of course Kcit] ri) = for t < n and for n slightly smaller than t it is clear that 
Kc{t;n) decays exponentially in n. We also know that (K c (t; 0)) ~ t" ds/2 . We 
believe that (Kc(t;n)) is a decreasing function of n for a fixed £ but a proof does 
not seem to be straightforward. 

7.3 First passage 

Let qc(t; n) be the probability that a walk which leaves the root at time t = on 
a comb C hits the vertex n on the spine for the first time at time t. We define the 
corresponding generating function as 

oo 

U c {x;n) = ^2{l-x) t/2 q{t;n). (169) 
t=o 

The mean first passage time at n is defined as the quantity 



°° d 
tc{ri) = ^2tq c {t;n) = -2 —U c {x;n) 

t=o 



(170) 



x=0 



One can calculate the probability <ZjVoo(^ n ) explicitly by elementary combinatorics 
with the result 

Un x (x; n) = — 1 (171) 
cosh(m oc (a;)n) 
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which implies that the mean first passage time on the half line is too (n) = n 2 . 

In this paper we shall not attempt to make a full calculation of the averaged 
probability distribution (qc(t;n)) but rather be content with showing that all the 
average mean first passage times are infinite on random combs with a < \. For this 
purpose it clearly suffices to show that (ic{2)) is infinite. The physical reason for 
this is easily seen to be that if the teeth in a random comb are sufficiently long to 
shift the spectral dimension away from 1 then the average time spent in a random 
tooth is infinite. 

Let us first assume that we have a random comb where there is a non-vanishing 
probability p that we have an infinite tooth at each vertex on the spine. Then 



where the lower bound is obtained by taking only into account combs with an infinite 
first tooth and restricting the attention to walks that wander into the first tooth and 
proceed directly to the vertex 2 when they return from the first tooth. Similarly, if 
the teeth are finite with probability distribution 



Using (fr2|) the right hand side in (|173|) is easily seen to be bounded from below by 



which is infinite if a = (70 + l)/4 < 1/2 where 70 is defined by (JED). 

8 Discussion 

Heat kernels on graphs and Riemannian manifolds have been extensively studied 
by mathematicians, see, e.g., (221123] and references therein. Much of this work is 
aimed at establishing the connection between pointwise behaviour of the heat kernel 
and geometrical properties of the graphs and manifolds. The most relevant results 
from our point of view are inequalities which, in our notation, are written 




(172) 




(173) 



00 





2d H 



<d s < d H , 



(175) 



l + d H 
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valid for graphs where the Hausdorff dimension exists and is finite, see [22] Theorems 
2.2 and 2.3. The Hausdorff and spectral dimensions calculated in this paper are 
readily seen to satisfy these bounds, some saturate the lower bound, others saturate 
the upper bound and some are in between, see Table 1. 

It is not understood in detail which properties of a graph cause the spectral and 
Hausdorff dimensions to differ. For the special and rather simple example of combs 
it is easy to see that if we have mostly short teeth, then the spectral and Hausdorff 
dimensions are both equal to 1. As the teeth grow, both dimensions grow but the 
Hausdorff dimension grows faster in general. It would be interesting to relate the 
two dimensions to the distribution of the orders of vertices. For more general graphs 
the connectivity clearly plays a role, not only the order distribution. 





d H 


d s 


2d H 
l+d H 


Random tooth spacing 


2 


3 
2 


4 
3 


Random tooth length a > 2 


I 


I 


I 


Random tooth length I < a < 2 


3 — a 


4-a 
2 


6-2a 
4-a 


Growing tooth spacing 


2+a 
1+0 


3+a 
2+a 


4+2a 
3+2a 


Growing tooth length a > 2 


2 


3 
2 


4 
3 


Growing tooth length 1 < a < 2 


2 


2(l+a) 
2+a 


4 
3 


Growing tooth length < a < 1 


I + a 


2(l+a) 
2+a 


2+2a 
2+a 


Random trees 


2 


4 
3 


4 
3 



Table 1. The spectral and Hausdorff dimensions discussed in this paper. 

There is no analytical understanding of the spectral dimension of random sur- 
faces and higher dimensional random manifolds. The spectral dimension of random 
surfaces is believed to be 2 jSJ |S] while the Hausdorff dimension is known to be 4 [21] , 
see also [23 [2E] • It has been shown recently [23 [2H] that the generic structure of 
infinite planar random surfaces (triangulations) is analogous to that of the random 
infinite trees discussed in Section 5. If we take such a surface S with a marked vertex 
and look at the boundary of a ball B, the boundary will have a number of disjoint 
components and with probability 1 only one of these components bounds an infinite 
subsurface of S. This means that we can view the infinite planar random surface 
as a tube with finite size outgrowths (baby universes) which are in fact distributed 
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in a simple way analogous to (J91j) . The tube and the outgrowths on the random 
surface correspond to the spine and the teeth of the random comb. Whether this 
picture allows us to obtain a rigorous control over the spectral dimension of random 
surfaces remains to be seen. 
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Appendix 1 

In this appendix we compute the first return generating functions for a number 
of simple combs. The method is essentially to solve the recursion relation (JJJ) in 
simple cases. 

We begin by evaluating the first return generating function Pe(x) for finite teeth. 
We note that P\ (x) — 1 — x and 

p ' w = 2^ra (176) 



Let us define 



A, = (177) 
It follows from (|176|) that satifies the recursion relation 

Writing 

A = 2 ~ P °° and X^ = A7 1 (179) 

-Too 

we see that ()178j) can be written 

X t = ApQ_! - 1) (180) 

which has the solution 

X l+l = A t X 1 -A^-^. (181) 
Inserting the values of A and X^ and doing some algebra leads to the desired result 

(H2D- 

Next, let us consider the comb C(L) which has infinite teeth at 1,2, ... ,L but 
no other teeth. Then, by (0), 

1 — x 

3 - Poo - Pc{l-i)(x) 

Defining 



P C (L)(x) = = = — . (182) 



^ = / C p = % (183) 

O -*oo 



we find from (11821) that Ex satisfies the recursion relation 



3 — -Poo — -P* 1 — -Ex-i 
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which is of the same form as ()178j) . Hence, by the same reasoning as that leading 
to (dHU) we find that 

E l = i r>L (185) 

BLE- 1 - B—^- 
1 — B 



where 



^ — P — P 

B = g — . (186) 



Pom - P*= , , Zvr.* 7v (187) 



Hence, 

P — P 
1 + F(5 L - 1 
with 

Noting that 5 = 1 + 2a; 1 / 4 + 0[\fx) as x — > we see that 

F = ^ + 0(a; 1/4 ). (189) 

The comb C(L) clearly has spectral dimension 1 since it only has a finite number of 
teeth. However, for our application in Section 6, we are interested in the behaviour 
of Pc(L) (%) as x — ► with L behaving like a negative power of x. 

Let us now assume that L = [x' 13 } with < (3 < |. It follows by expanding out 
the denominator in (J187)) that 

P c(L) (x) = 1 - to*-? + o(x*~ p ). (190) 

We conclude that 

l-P c(L )(x) ~ x h-P (191) 

which was needed for the inequality (|12Hj) . 

The final comb we consider in this section has no teeth at vertices 0, 1, ... k — 1 
but teeth of length I at k, k + 1, . . .. Denote this comb by C(k,£). We will be 
interested in the limit of the first return generating function on C(k,£), denoted 
Pk,e, when k and £ tend to infinity as a; — > 0. The logic of the calculation is very 
much the same as above. 

Let us denote P\i by P# and note that this function satisfies the recursion 
relation 
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k 



Figure 4: The comb C(k,£) with k = 5 and i = 4. 



which is easily solved with the result 
With A defined as in (|179j) one finds 



i - p e ( x ) + -(l - p e ( x )y + x . 



(193) 



? — P 

* 1 CO 

Pk,£ — Poo 

cf. pSljl . This can be written 

Pki — Poo 



A 



? — P 

-1 * 1 oc 
P# - Poo 



-A 



I- A 



k-l 



1 - A 



(194) 



2(P # -P 0O )(l-P 0I 



2(1 - Poo) + (2 - Poo - P#){A k ^ - 1) ' 

We now wish to find the asymptotic behaviour of Pk,i{x) as x — > with fc 
and £ = [A; a ]. Assume that < /3 < ~ and recall that 

1 + v^ 



(195) 



;r 



.4 



1- 



It follows that 



A fc - 1 = 2 + 0(2^) 

Note that (|12|) can be written 

v4* — 1 
1-P(x) = v^- 



+ r 

so if /3a < I and we use (I197j) with /3 replaced by /3a we find 

1 -P e ( x ) = x 1 -* 3 " + o(x 1_/3 °). 
In the case /3a > ~ it is not hard to see that 



(196) 



(197) 



(198) 



(199) 



1-P e (x) = ^x + E(x) 



(200) 
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In the limiting case (3a = \ an easy calculation yields 



2 

e 2 - I 

e 2 + 1 



l-P^(a;) = :T - rT V^ + o(V^)- (201) 



The next step is to find how P#{x) behaves as x — > and then we can infer the 
behaviour of Pk/{x) from ()195j) . It is convenient to split the argument into two 
cases. 

Case 1: (3a>\ 

It is clear from (|193|) that P#(x) = 1 — cx? + 0(2:4). Using this we find that 

P k/ {x) - P M (x) = - C f + ' + ? ' (202) 

If (3 < ~ then we find 

PkA x ) = P °o(x) - ex* + o{xi) (203) 

and if (3 > \ then 

P k Ax)=Poo{x)-x> 5 + o{x p ). (204) 

In the crossover case (3 = \ we find ()203|) with the constant c replaced by c(l + c) _1 . 
We remark that this calculation is insensitive to the value of a. 

Case 2: (3a < \ 

Using (USSD and (M we obtain 

, ^ l-0a , 1 /3 d 

p # ( x ) = + (a;— ). (205) 

It follows that 

2-fia 2-)3a 

P M (x) - Poo(x) = - * \ + (206) 

V / X + X 1_ ^ _ 2 /3a (l + O(l)) 

There are again two cases to consider. If 

\<l-(3-\(3a (207) 

then 

^(a:) = P oo(x) - ca^ 1 "^ + oix^ 1 -^) (208) 
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and c = \ if equality holds in ()207|) but c = 1 otherwise. If | > 1 — (3 — ^f3a then 

PkAx)=Poo(x)-x l3 + o(x> 3 ), (209) 



i.e., 

^/(z) = Poo(x) - ca; 5 + o(x s ) (210) 
where 5 = max{i(l — (3a), [3}. 

Appendix 2 

In this appendix we show that for Lyfx small 

i 

P #L (x) = l--^ + 0(v^), (211) 

and hence, for L = with (3 < |, we obtain 

P^ (x) = 1 - + o(x2 + ^) (212) 

which is the result needed in Sections 3 and 6. 

Let Rl denote the generating function for first return to the root for walks that 
do not move beyond the (L — l)st vertex on the spine, i.e., they do not reach the 
vertex where the first tooth appears. Let denote the two-point function defined 
by the sum over all walks from the root to L which do not return to the root and 
stop the first time they meet L, i.e., Tl(x) = Gq(x;L) where C = N^. Then by 
decomposing the walks that contribute to P* L we obtain the equation 

r 2 

P. =R L + £ . (213) 

° r oo r *L n L 

It is straightforward to solve this to obtain 
3-P 



oo 



± o V ( 3 - p oo) 2 - 4(r| + P L (3 - Poo - R L )). (214) 



^ 2 2 

Since P* L (x) < 1 we must choose the — sign in (|214|l . The calculation of Rl is 
similar to that of Pi in Appendix 1 and gives 

Kl W-^ x} (i + ^-(i-v^) L [ } 

so in particular 

Rl{x) = ^-^ + 0(xL). (216) 

Lj 
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Using an argument similar to the one leading to (|2"ij) we see that 

T L (x) = (1 - x)-^ L+1 R L (x)R L ^(x) . . . R 2 (x) (217) 

for L > 2. If L = 1 then the product of the R factors in ()217|) is equal to 1 by 
definition. Rearranging we find 

T L (x) = (1 - xfl \ ^ -=- (218) 

(1 + y/x) L - (1 - V X ) 

SO 

r L = ^ + 0(Lx). (219) 
Noting that the expression under the square root in ()214|) can be written 



(3 - Poo - 2R L + 2r L )(3 - Poo - 2P L - 2T L ), (220) 
we obtain (JTTT]) by inserting ffTTl) . lj2T7)|) and (jTTTIjl into (gUD - 
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